Higher harmonics increase LISA's mass reach for supermassive black holes 
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Current expectations on the signal to noise ratios and masses of supermassive black holes which 
the Laser Interferometer Space Antenna (LISA) can observe are based on using in matched filtering 
only the dominant harmonic of the inspiral waveform at twice the orbital frequency. Other harmonics 
will affect the signal-to-noise ratio of systems currently believed to be observable by LISA. More 
significantly, inclusion of other harmonics in our matched filters would mean that more massive 
systems that were previously thought to be not visible in LISA should be detectable with reasonable 
SNRs. Our estimates show that we should be able to significantly increase the mass reach of 
LISA and observe the more commonly occurring supermassive black holes of masses ~ 10 s Mq. 
More specifically, with the inclusion of all known harmonics LISA will be able to observe even 
supermassive black hole coalescences with total mass ~ 1O 8 M0(1O 9 Mq) (and mass-ratio 0.1) for a 
low frequency cut-off of 10 _4 Hz (10 _5 Hz) with an SNR up to ~ 60 (~ 30) at a distance of 3 Gpc. 
This is important from the astrophysical viewpoint since observational evidence for the existence of 
black holes in this mass range is quite strong and binaries containing such supermassive black holes 
will be inaccessible to LISA if one uses as detection templates only the dominant harmonic. 

PACS numbers: 04.30.Db, 04.25.Nx, 04.80.Nn, 95.55.Ym 



I. INTRODUCTION 



A. Supermassive black hole binaries and LISA 



There is strong observational evidence for the existence of supermassive black holes (SMBHs) with masses in the 
range of 10 6 Mq-10 9 M© (see e.g. Ref. [l[ and references therein) in most galactic nuclei @]. Therefore, mergers of 
galaxies, as evidenced by high-redshift surveys, should give rise to binaries containing SMBHs. Late stage evolution of 
a SMBH binary is dictated by the emission of gravitational radiation. The resulting loss of orbital energy and angular 
momentum would lead to the coalescence of the two holes. Indeed, X-ray observations have revealed the existence of 
at least one such system that would coalesce within the Hubble time Gravitational waves (GW) emitted in the 
process could be detected by the planned Laser Interferometer Space Antenna (LISA) Q. 

Observation of SMBH binaries at high redshifts is one of the major science goals of LISA. These observations will 
allow us to probe the evolution of SMBHs and structure formation Q and provide an unique opportunity to test 
General Relativity (and its alternatives)in the strong field regime of the theory @, 0, 0, S, lld l. Observing SMBH 
coalescences with high (100-1000) SNR §1, [t| is crucial for performing all the aforementioned tests. 



B. Restricted Vs Full Waveforms as Search templates in LISA 

Motivated by the fact that matched filtering is more sensitive to the phase of the signal than its amplitude [Til ] . 
search algorithms so far have deployed a waveform model involving only the dominant harmonic (at twice the orbital 
frequency), although the phase evolution itself is included to the maximum available post-Newtonian (PN) order 
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(currently 3.5PN, for non-spinning systems [Hi EH)- Waveforms in which all amplitude corrections arc neglected, but 
the phase is treated to the maximum available order, are called restricted waveforms (RWF) and these are what are 
used so far in the analysis of data from ground-based detectors [3, EEI, EH, E3 ■ This paper will consider the advantage 
of using the full wave forms (FWF) in the context of LISA. 

LISA is designed to detect gravitational waves in the frequency-band 0.1-100 mHz. This frequency range determines 
the range of masses accessible to LISA because the inspiral signal would end when the system's orbital frequency 
reaches the mass-dependent last stable orbit (LSO). In the test-mass approximation, the angular velocity ulso a t 
LSO is given by wlso = 6~ 3//2 Af _1 , where M is the total mass of the binary. Search templates that contain only 
the dominant harmonic cannot extract power in the signal beyond / LS o = wlsoA - 4.39(Af/lO 6 M )- 1 mHz. This 
further implies that the frequency range [0.1, 100] mHz corresponds to the range ~ 4.39 x [10 4 , 10 7 ]M Q for the total 
mass of binary black holes that would be accessible to LISA 1 . However, as Table 1 of Ref. [l[ would reveal, there is 
observational evidence for the existence of many SMBHs whose masses are of the order of 10 8 -10 9 M Q . LISA will be 
unable to observe binaries containing SMBHs in this mass range if it used as search templates waveforms containing 
only the dominant harmonic. 

Inclusion of higher-order amplitude terms in the waveform introduces the following two new features: (i) appearance 
of higher harmonics of the orbital phase and (ii) PN amplitude corrections to the leading as well as higher harmonics of 
the orbital frequency For example, at 0.5PN order, which is the first-order correction, there are two new harmonics VP 
and 34', where *f> is related to the orbital phase of the binary as in Refs. [H, EH- More interestingly, in the expressions 
for the 'plus' and 'cross' polarizations, all odd harmonics of the orbital frequency are proportional to where Sm is 
the difference in the masses of the binary components (see Eq. (5.7)-(5.10) of Refs. jlj|). Another important feature 
of the full waveform is that the (2n + 2) th harmonic first appears at the n th PN order in amplitude 2 . For example, 
the fourth harmonic first appears at 1PN, and has PN amplitude corrections to its dominant term at 2PN and 2.5PN 
(see Refs. [H,E1 for details). 

Early investigations on the importance of amplitude-corrections to search templates were carried out by Sintes and 
Vecchio [13,[2l|. Their study used only the first-order correction at 0.5 PN order. They concluded that the addition of 
the amplitude terms in the waveform did not improve the accuracy in the estimation of source's angular position and 
the distance, whereas the estimation of the chirp and reduced masses could be 10 times better when compared to the 
RWF. Recently, in the context of ground-based detectors, Van Den Broeck and Sengupta [13, HH, [24| examined the 
implications of going beyond the restricted PN approximation and employing instead the full waveform [Tel Il9| . The 
two main implications of the comprehensive analysis in Refs. [22j, |23j for terrestrial GW detectors may be summarized 
as follows: 

1. For binary neutron stars and stellar mass black holes, restricted waveforms over-estimate the SNR as compared 
the full waveform. 

2. The use of the full waveforms significantly increases the mass-reach of second and third generation detectors, 
advanced LIGO and EGO being able to observe systems with total mass ~ 400M Q and a third generation 
detector as high as 10 3 M©. 

In the present paper, we study in the context of LISA the implication of using templates based on the FWF (i.e. 
including all known harmonics of the orbital phase and all known amplitude corrections in the GW polarisations). 
Coalescences of SMBH binaries with masses ~ 10 8 ~ 9 M© will not be observable by LISA if one uses only templates 
based on the RWF. Using templates based on amplitude corrected full waveforms, instead of the usual restricted 
waveforms, will enable LISA to observe coalescences of SMBH binaries with total mass ~ 10 8 Af Q (10 9 Mq) if the 
lower frequency cut-off LISA can achieve is ~ 10 _4 Hz (10 _5 Hz). 

The rest of this paper is organized as follows: In Section[Hj we give the FWF in the frequency domain, by taking into 
account the orbital motion of LISA around the sun and its changing orientation. Section IIIII discusses the results of 
our investigations where we compare the performances of the amplitude-corrected waveforms at different PN orders in 
terms of their mass- reach and distance- reach and correlate it to the 'observed' spectrum in LISA. Section fVl concludes 
with a brief summary of the main results and assumptions underlying their derivation. 



1 Although, binaries lighter than 1O 4 M0 would, in principle, evolve through the LISA band they would not be luminous enough to be 
visible in LISA unless they are close-by. 

2 The 0.5PN term is an exception to this and also introduces a harmonic at the orbital frequency apart from the one at thrice the orbital 
frequency. 
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II. TEMPLATE WAVEFORMS FOR LISA 
A. Amplitude corrected waveform 

For non-spinning binaries in quasi-circular orbits inspiralling due to radiation-reaction, waveforms were computed in 
Refs. [IH, [lj| up to 2.5PN order in amplitude and 3.5PN in phase [IH[l3||. This waveform h(t) is a linear combination 
of sine and cosine functions of multiples of the orbital phase ^(i). The expression for the 2.5PN polarization contains 
the first seven harmonics of the orbital phase, the dominant harmonic being the one at twice the orbital phase. The 
signal depends on the following parameters: the luminosity distance to the binary, m the total (red-shifted) mass, 
v the symmetric mass-ratio (reduced mass divided by total mass), the spherical polar angles (#,</>) determining the 
direction of the "line-of-sight" , the inclination angle l of the angular momentum L of the binary with respect to the 
direction opposite to the line-of-sight, and the polarization angle ip which determines the orientation of the projection 
of L in the plane normal to the line-of-sight. 

We rewrite the waveform in terms of only cosines in a form similar to [23[ : 

2Mv 7 ^ 

h(t) = Y C Os[k^(t) + 4>{k, n /2)] ^ + \t), (2.1) 

L fc=ln=0 

where the coefficients A^ n /2) and <p(k,n/2) are functions of (v,9,(l),'ip,i), and x(t) = (27rAfF(i)) 2 / 3 is the post- 
Newtonian parameter with F(t) the instantaneous orbital frequency. Terms x n ^ 2+1 (t) A^ k n / 2 ) and 4>{k,n/2) are 
the wave amplitude and polarization phase, respectively, corresponding to the k th harmonic and (n/2) th PN order. 
We call the coefficients A^. n /2) the polarization amplitudes. The orbital phase ^(t) is a PN series in x, which, in 
the case of non-spinning binaries, is known to 3.5PN order [l2l. [T3|. For a non-spinning source and a detector whose 
position and orientation are almost constant during the time of observation of the signal, all the above mentioned 
angles are constants. For ground-based GW detectors dealt with in Ref. @, one is in this situation. 

B. Amplitude corrected waveform including modulations due to LISA's orbital motion — Time Domain 

LISA will be able to observe many sources from their early stages of inspiral and most would last for a pretty long 
time. We shall only consider binary sources that last for a year or less before merger. Since the LISA plane is tilted 
by 60° with respect to the plane of the ecliptic, during the course of its heliocentric orbit its orientation and position 
varies periodically, with a period of one year and the signal in Eq. (|2.1[) will suffer additional amplitude and phase 
modulations. Thus in the case of LISA the angles 9, </>, and ip (but not i) appearing in Eq. 12. H are functions of time. 
To proceed further, in the frame of a non-rotating observer fixed to the solar-system barycentre, we denote by the 
location of the source on the sky by the spherical polar angles 9s and <ps and the orientation of the source by the 
spherical polar angles 9l and <pL determining the direction of the orbital angular momentum L of the binary. The 
transformation between the fixed set of angles 3 (6s, <t>s, 9l, <Pl) and the time-dependent angular coordinates of the 
source (9(t), 4>(t), tp(t),i) as measured by LISA are given in Ref. (2f|. 

Generalizing, Ref. [25[ from the RWF to the FWF, the signal as seen in LISA is of the form, 

h(t) = — -p-^ Y A (k,n/2)(t) cos[fc*(t) + (M/a) (t) + k</> D (t)] x% +1 (t) . (2.2) 

The PN parameter x(t) appearing in Eq. (|2.2p is still equal to (2irMF(t)) 2 / 3 , where F(t), however, is the orbital 
frequency as measured by a non-rotating observer located at the solar-system barycentre. The term (f>D{t) is the 
Doppler phase [25[ , accounting for the phase difference of the gravitational wave- front between LISA and the solar- 
system barycentre. The time-dependence of 4>D(t) is due to the orbital motion of LISA about the barycentre. It is 
given by 

t/> D (jt) = 2irF(t)Rwi6scoB[4>(t)-tl>s], (2.3) 

where R = 1 AU and 4>(t) is the angular position of LISA with respect to the barycentre given by <f)(t) = 2tt ^, T 
being equal to one year. 



This is a different notation from l2- r | . where the source angles measured in the fixed barycentre frame are denoted by (0$, 0g, 9^, c/>l) 
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C. Amplitude corrected waveform including modulations due to LISA's orbit — Frequency Domain 

The above waveform is valid in the adiabatic regime, where the radiation-reaction time-scale is much larger than 
the orbital time-scale. We also note that the additional amplitude and Dopplcr modulations in the waveform for LISA 
vary on time-scales of 1 yr (i.e. ~ 3 x 10 7 s), while LISA can observe orbital periods at most up to 2 x 10 5 s, (i.e. 
gravitational wave frequencies of order 10 Hz.). Consequently, the Doppler modulations change much more slowly 
(a hundredth) than the orbital phase. This permits the use of the stationary phase approximation (SPA) to obtain 
an analytical form for the Fourier transform (FT) h(f) of the signal: 

V32Mv * * A (fc , Tt/2) (t(//fc))xf+ 1 (t(//fc))e-^^/ 2 )W/A)) 

Kf) - -7T-FT 1^ -* lJ - LJ j exp[zV/(t(//fc))] , (2.4) 

& fe =i n=0 2^JkF(t(f/k)) 

where an over dot denotes derivative with respect to time and tpf(t(f/k)) is given by 

MKf/k)) = 27r/i(//fc) - fc*(f (//*)) ~kcp D {t{f/k)) - tt/4. (2.5) 

The PN expansions for t(F), ^f(F), F(F) are given in Ref. [l2j . The expression in Eq. ()2.4|) within the summation over 
k represents the FT due to the k th harmonic. It should be noted that the term F may be treated in different ways 
that could lead to numerically different results. In a numerical treatment, for instance, one could avoid performing a 
further re-expansion. Alternatively, one could re-expand the denominator in the amplitude and truncate the resulting 
expression at the n th PN order, to obtain the nPN amplitude-corrected waveform. Ref. [22|, [23| choose the latter and 
we follow them in this work. 

Radiation reaction results in an increase in the orbital frequency F(t) which will ultimately drive the system beyond 
the adiabatic inspiral phase and the inspiral waveform given above will no longer be valid. In the first approximation 
this is expected to occur when the orbital frequency F(t) reaches FLso ~ the orbital frequency of the LSO of a 
Schwarzschild solution with the same mass as the binary's total mass M, 

F LSO = (2 7r6i M)- 1 . (2.6) 

Thus, we truncate the signal in the time domain at a time t LSO ; given implicitly by F(^lso) = Flso- In the SPA, 
the main contribution to the FT of the k th harmonic at a given Fourier frequency /, comes from the neighbourhood 
of the time when the instantaneous value of the k th harmonic sweeps past /. Thus the fc th harmonic is not expected 
to contribute significant power to the FT for frequencies above AiFlsOi if the signal is truncated in the time domain 
beyond iLSO- This motivates the truncation of the FT due to the k th harmonic at frequencies above k -FLso by a step 
function 9(kFj J so ~ f) [6( x ) = L for x>0, and for negative x\. 

III. SIGNAL TO NOISE RATIOS IN LISA WITH HIGHER HARMONICS 

In this Section we investigate the effect of the higher harmonics in LISA observations of supcrmassivc black hole 
binaries. The LISA waveform discussed in the previous Section will be used for the analysis. 

Given a waveform h, the best signal-to-noise ratio (SNR) achieved using an optimal filter is given by p[h] = (hlh) 1 / 2 , 
where ( . | . ) is the usual inner product in terms of the one-sided noise power spectral density Sh{f) of the detector. 
With the convention for Fourier transforms, x(f) = x{t) exp(— 2irift) dt, the inner product is given by: 

{Av)si j'-^fMl d , (3.!) 

For an optimal filter, which maximises the overlap of the signal with template, one can write 

We use the non-sky-averaged noise-spectral-density as given in Eqs. (2.28)-(2.32) of Ref. 0. 
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waveforms, for two archetypal binaries, (10 s , 10 )M@ (left) and 2 x (10 , 10 )Mq (right). The sources are assumed to be at 3 
Gpc and their orientation with respect to the solar-system barycentre is chosen to be 8s = cos -1 (— 0.6), <f>s = 1, 0l = cos -1 (0.2), 
4>l = 3. The spectrum is much more complicated and highly oscillatory for the FWF than for the RWF, because of interference 
between various harmonics. The higher frequency reach of the FWF is due to presence of higher harmonics as apparent in the 
figure. The spectrum for the system in the left panel sharply rises at a frequency ~ 2 x 10 _3 Hz. Beyond this frequency, the 
effective LISA noise decreases sharply with increasing frequency (as there are fewer galactic binaries per frequency bin) leading 
to the observed increase in the spectrum. 



A. Choice of frequency cutoffs f en d,fs 

The upper limit of integration / en( j is taken to be the minimum of 7 -Flso and 1 Hz, the latter being a conventional 
upper cut-off for the LISA noise curve. The lower limit / s is chosen assuming LISA observes the inspiral for a 
duration Ai b s before it reaches the LSO. Using the quadrupole formula, we find that the orbital frequency at the 
epoch t LSO - Ai obs is given by 



F(t 



LSO 



Aiobs) — 



S< ) 



(1 



256 v 



5M ^obs v LSoJ 



(3.3) 



where i>lso is the orbital velocity and ^lso the epoch at which the orbital frequency reaches the value -FLso- We 
take vlso = l/\/6, the orbital velocity at the LSO in the case of a test mass orbiting a Schwarzschild black hole. We 
designate Ffaso — Ai b s ) as Fj n . Thus the fc th harmonic will have a frequency fc-F] n , Ai b s before £lso- The above 
formula reduces to the simpler form in Ref. [|| as wlso ~~ * °°- F° r the mass values explored in this work there is no 
significant dependence of the results on this choice. In all our calculations we take At b s to be one year. 

The lower cut-off for the fc th harmonic should be the maximum of the lower cut-off of LISA (10 -4 Hz) and k F ln 
and simply implemented by truncating the waveform due to the fcth harmonic by another step-function 6{f — k Fi n ) 
and choosing / s to be 10 _4 Hz. It is worth noting that the fc th harmonic probes a larger interval of the frequency 
domain i.e. fc(i 7 Lso — F[ n ) relative to the fundamental harmonic. For brevity, we refer to this as the span of the fc th 
harmonic. 

There is a caveat with regard to the use of higher harmonics that is worth mentioning: In the time-domain the 
waveform should begin when the highest harmonic reaches the lower cutoff. This has an implication on data analysis 
as the templates will be an order-of-magnitude longer than before. Thus, it might be sensible to use higher harmonics 
only in the case of higher masses. 



B. Observed signal spectrum with LISA 



To get some insight into the effect of higher harmonics via amplitude corrections let us first look at the SNR 
integrand, i.e., the "noise- weighted signal power" per unit logarithmic frequency interval pfjj . Rewriting the expression 
for the SNR as 



f\Hf)\ : 
s h (f) 



dln(f), 



(3.4) 
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TABLE I: SNRs due to successive PN amplitude-corrected waveforms, with phase corrections to 3.5 PN order in all cases. The 
orientation of the source with respect to the solar-system barycentre is chosen to be 9$ = cos _1 (— 0.6), <f)s = 1, 9l = cos (0.2), 
4>l = 3. For the (10 6 — lO 7 )Af0 binary system, all harmonics enter deep into the sensitive part of the LISA bandwidth. Apart 
from an increase at 0.5PN, we see a consistent reduction in the SNR on inclusion of higher PN order amplitude corrections. For 
the (5.5 x 10 6 , 5.5 x 1O 7 )M0 binary system, the second harmonic fails to enter the LISA bandwidth, while the third harmonic 
spans a small insensitive region. Thus the SNR due to the RWF is zero, while the SNR due to the 0.5PN waveform is smaller 
than the SNRs due to higher order PN terms. Both sources are at a distance of 3 Gpc. 



PN 
order 


SNR 


(10 6 - 10 7 )A/ Q 


5.5 x (10 6 - 10 7 )M© 





924.48 





0.5 


1025.8 


211.98 


1 


928.48 


343.17 


1.5 


869.78 


319.34 


2 


824.65 


266.65 


2.5 


809.51 


277.34 



the quantity of our interest is 

which is designated as the "observed spectrum" , following |24j . The observed spectrum is plotted versus frequency for 
given masses in Fig. [TJ As is the case for ground-based detectors (24[, the spectrum due to the FWF has a lot more 
structure and is highly oscillatory because of interference between various harmonics. For the (10 5 , 10 6 )Af Q system, 
the mass being low, the second harmonic and hence the RWF extends up to frequencies ~ 2 x 10 -3 Hz, where LISA 
is most sensitive. This leads to a rapid increase in the observed spectrum in this frequency region. The spectrum due 
to the FWF, containing higher harmonics continue beyond the RWF into the most sensitive part of the LISA band. 
For the 2(10 6 , 10 7 )M Q system, the frequency span of the second harmonic is small and the sensitive region of the 
LISA band lies beyond its maximum reach. 



IV. THE EFFECT OF HIGHER HARMONICS 



Following the analysis of Ref. [22], |23[, we classify the sources into two types: In the first category are sources for 
which the dominant (second) harmonic has a large frequency span in the LISA band. The second category on the other 
hand comprise sources whose dominant harmonic fails to enter the LISA bandwidth but the higher harmonics do. 
Since the upper cut-off frequency for each harmonic is inversely proportional to the total mass (from the expression 
for FLso), we note that the sources of the first type will have total mass less than some value which we call the RWF 
mass-reach, the maximum mass detectable by the RWF, while the second type will have masses greater than this 
value. The condition that the upper cut-off of the dominant harmonic is less than or equal to the lower cut-off of 
LISA (i.e., by the inequality 2Flso < fs) determines the RWF mass-reach. The choice of / s for the LISA mission is 
still not clear and theoretical implications of this choice are explored in e.g. Ref. (27} . For / s in the range [10 -5 , 10~ 4 ] 
Hz the RWF mass- reach varies over the range [4.39 ,43.9] x 10 7 M Q , the lower end of the mass range corresponding 
to the higher end of the frequency range. 



A. How higher harmonics affect signal visibility 



In Fig. [2] we plot the SNRs computed using the restricted (RWF) and full (FWF) waveforms as a function of 
the binary's total mass for two values of the mass ratio 4 . We first consider systems whose total mass is less than 



4 Our codes are calibrated by reproducing the results of |22|. l23l . which considers ground-based detectors, and also of 28], which computes 
SNRs in LISA using RWF. 
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FIG. 2: SNR versus total mass for mass ratios of 0.1 and 0.01. The figure shows that apart from the dips due to white-dwarf 
confusion noise, for mass values where the RWF enters the LISA band, the corresponding SNR is consistently more than 
the SNR produced by the FWF. However, for mass values where the second harmonic terminates before it reaches the LISA 
bandwidth, the FWF which has higher harmonics that enter the LISA band produces significant SNRs. The frequency reach 
of a harmonic depends only on the total mass, and thus the mass reach of the FWF is independent of the mass-ratio. For 
more asymmetric systems, the magnitude of the SNR is low for all masses both for the RWF and the FWF. Sources are at a 
luminosity distance of 3 Gpc with fixed angles given by 6s = cos _1 (— 0.6), cf>s = 1, 6l = cos _1 (0.2), 4>l = 3. 

4 x 1O 7 M . For these systems, the SNRs computed using the two different approximations agree with each other 
to within 10%, with the RWF over-estimating the SNR, when compared to the FWF. in most of the range. This is 
explicitly shown for a (10 6 , 10 7 )M Q binary in the first column of Table U Indeed, but for the slight increase in SNR 
as we go from 0PN to 0.5PN, we find a steady decrease as one increases the PN order of the amplitude correction. 

The reduction in SNR at higher PN orders can be understood by studying the structure of \h(f)\ 2 , the numerator 
in the integrand of the SNR in Eq. (|3.2[) . There are basically three types of terms: 

1 . direct terms in which the phases in Eq. (|2.4[) cancel 

4*,n/2)W/A))/" 5 

2. interference terms between different PN corrections of the same harmonic, 

A (k ,m/2)(t(f /*)) A (Kn/2) {t{f ■/*)) /-5 (Mf)^ cos[0 (W2) W//fc)) - (fc ,„/2)W// fc ))] 

3. harmonic mixtures 5 which are terms containing the interference between different PN corrections of different 
harmonics, e.g. the m/2 th PN correction of the fc th harmonic and n/2 th PN correction of the I th harmonic. 

A {k>mm (Af/k)) A M2) (t(f/l)) f-i (M/) 2 * 2 ca#/(t(//fc)) - (W2) (t(f/k)) - rf> f (t(f/l)) + (i ,„ /2) (t(f/l))\ 

where ipf(t(f/k)) is given by Eq. (|2.5p . 

All these terms are scaled by where M. = M v 3 / 5 is the chirp-mass. (Additional multiplicative factors have 

been omitted in the above expressions, among which are the step-functions mentioned earlier and PN expansion 
coefficients of the denominator of the Fourier amplitude in Eq. (|2.4p . the latter being time- independent.) 

B. The effect of higher harmonics in ground-based detectors 

Before we explain the SNR trends in the context of LISA, we mention that for ground-based detectors a similar 
effect was found in Ref. [23| for a different but corresponding mass region. The lower cut-off for a typical ground-based 



5 We use the term 'harmonic mixtures' at the risk of being mistaken to the well-known 'harmonic mixing' in music. Our use of the phrase 
'harmonic mixtures' is simply to convey the physical effect of the interference between different harmonics 
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FIG. 3: Variation of polarisation amplitude of the RWF with frequency and time (inset). The inset, plotted over a duration of 
two years clearly shows periodicity due to LISA's orbital motion around the Sun. The binary mass, (10 6 - lO 4 )Af , has been 
chosen such that it can, in principle, be observed for two years. The plot in the frequency domain shows that the variation 
of the polarisation amplitude is confined to a very small part of the frequency span of the dominant harmonic, and essentially 
behaves as a constant in the frequency domain. 



detector, say Advanced LIGO is 20Hz, and the effect of higher harmonics is seen for masses less than ~ 220M Q . In 
that case, as mentioned earlier, the polarisation amplitudes and phases are constants. The RWF contains only the 
Newtonian term of the second harmonic and thus \h(f)\ 2 consists of a single direct term with n = and k = 2. 

With the inclusion of higher-order amplitude terms in the waveform, PN corrections to the dominant harmonic, and 
higher harmonics and their PN corrections, also contribute to the SNR. In other words, the signal power spectrum 
IM/)| 2 wm contain all three types of terms discussed before. From the form of the direct terms, it is evident that 
their contribution to the SNR will be positive definite. We also note that, for ground-based detectors, the frequency 
dependence of the direct and interference terms will just be a power law. However, the sign of the interference 
terms (and consequently their contribution to the SNR) depends on the difference between the polarisation phases of 
different PN corrections for the same harmonic. Van Den Broeck and Scngupta showed that for a given harmonic, 
for all allowed values of the parameters (v, 9, (j>, ip, i), each PN correction is almost "out of phase" with both the PN 
correction preceding and succeeding it 6 . The resulting negative terms (representing destructive interferences) reduce 
the SNR as one includes higher PN amplitude corrections in the waveform. 

The third type of terms, harmonic mixtures, however, are highly oscillatory functions of the frequency, as the phase 
difference il>f(t(f/k)) — ipf{t(f/l)) between the k th and the I th harmonic become even or odd multiples of ir. As one 
integrates over /, these oscillations tend to cancel out, and thus the contribution to the SNR from these terms are 
numerically much smaller relative to the first two types of terms. 



C. Effect of higher harmonics for binaries with M < 4 x 10 t Mq 

In the case of LISA, because of the polarisation factors, the amplitudes of none of the three types of terms is a 
simple power-law in /. The periodic variation of, for example, -A( 2 .o) (period being one year) appears as an amplitude 
modulation ^4(2, o) W//2)) in the Fourier transform, where the argument t(f/2) of ^.2,0 is given by 

t(//2) = " 2567^5/3 7573 + PNcorrcctions - 

Hence, in the frequency domain A^.o) wn l undergo one complete oscillation as / varies from 2F m (see Eq. (|3.3jl ) to 
2-Plso- However, because of the inverse power-law dependence on /, the oscillation of -A( 2 .o) is confined to a small 
frequency interval above F[ n and remains fairly constant over a major portion of the frequency span 2(^30 — -Pin) 
(see Fig. [3]). For masses higher than the one shown in Fig. [3J this region of significant variation moves to the left of 
the figure. On including in our analysis the effect of detector sensitivity (weighting down by Sh(f)) this variation of 



Note, however, that Ref. [23ll . argues this in a somewhat different form. 



-• Newtonian 
•- 0.5PN 

1PN 
•- 1.5PN 

2PN 
-• 2.5PN 



Total Mass (in Msun) 



Newtonian 
0.5PN 
1PN 
1.5PN 
2PN 
2.5PN 



Total Mass (in Msun) 



FIG. 4: SNR versus total mass for successive PN amplitude-corrected waveforms and 3.5PN phasing. The left panel corresponds 
to a mass-ratio of 0.1 while the right panel plots the same quantities for mass-ratio of 1 (equal mass systems). The (2n + 2) th 
harmonic first appears at the nth PN order. For a binary of given total mass, the upper cut-off of the fc th -harmonic of the orbital 
frequency in the frequency domain is proportional to k and inversely proportional to the total mass. As the mass increases the 
upper cut-off for the 2 nd harmonic falls below the lower cut-off of the LISA detector, leading to a zero value of SNR due to the 
RWF. The higher harmonics still enter the sensitive bandwidth of LISA and higher PN order waveforms produce significant 
SNR. The 2.5PN waveform has the highest mass-reach, being 3.5 times the mass-reach of the RWF. In the equal mass case 
displayed in the right panel, the differences in harmonic content of different PN order waveforms are more pronounced, as odd 
harmonics are absent. Sources are at a luminosity distance of 3 Gpc with fixed angles given by 8s = cos~ 1 (— 0.6), 4>s = 1, 
6 L = cos" 1 (0.2), <j> L =3. 



A(2,o) gets damped out when one evaluates the integral in Eq. (|3.2[) . For masses satisfying 2F in <§C 10~ 4 Hz, the lower 
cut-off for LISA, this region of variation will fall below the LISA band. 

The polarisation phases determining the sign of the interference terms between the same harmonics also vary with 
/. However, as mentioned earlier, the phase relationships of the polarisation phases are independent of the parameter 
values. Thus the modulations which change the values of (6, <\>, ip, i) do not affect the trend of reduction of SNR with 
amplitude corrections. The Doppler modulations, which appear in only harmonic mixtures, are also not important as 
far as SNR is concerned. 

Finally, we would like to note an important point not explicitly mentioned in Ref . [23| . As the difference between the 
polarisation phases of successive PN corrections of the same harmonic tend to be nearly ir, alternate PN corrections 
necessarily interfere constructively. Hence there are positive contributions also from the interference terms. Now, the 
numerical value of the contribution to SNR from each of these terms depends on the magnitude of the polarisation 
amplitude and the power of (Mf). It can be checked that for all allowed values of (v,9,4>,ip,i) > the polarisation 
amplitudes are roughly of the same order in magnitude. Consider the Newtonian term of the dominant harmonic 
and its interference with the first two corrections at 1PN and 1.5PN. The Newtonian term will be out of phase with 
the 1PN term, but in phase with the 1.5PN one. The two corresponding interference terms will contain powers of 
(A//) 2 / 3 and {Mf), respectively, and since they have the same frequency span, the absolute numerical value of the 
contribution to SNR from the former will be more since (2ttM f) 1 / 3 will always be less than I/a/6. Numerical values 
of contributions from interference between higher PN corrections of the second harmonic successively decrease. The 
same argument applies for all the other harmonics, and thus, inclusion of amplitude corrections will lead to an overall 
reduction in SNR. 

The first column of Table |I] clearly demonstrates the effect of higher harmonics on SNRs. The increase in SNR 
for the 0.5PN waveform (with respect the RWF SNR) is also easily explained by noting that the 0.5PN correction 
only adds (apart from harmonic mixtures) two direct terms to \h(f)\ 2 , corresponding to the first and third harmonics 
(n = 1, n = 3). Clearly, from the discussion in the previous subsection, the 0.5PN waveform will have a higher SNR 
than the RWF, independent of the binary parameters. 

For 10 3 M ^ 10 5 M Q , the difference between the RWF and the FWF is not visible on the scale of Fig. [2] because 
for this mass range all the direct and interference terms corresponding to harmonics higher than the dominant ones, 
which are scaled by higher powers of (Mf), are negligible. 
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D. Visibility of systems with M > 4 x 10 7 M Q 

In their analysis of the implications of the FWF for ground-based detectors Van Den Broeck and Sengupta [22|, HH 
pointed out an interesting effect due to higher harmonics. An analogous effect is found in the case of LISA in spite 
of the additional amplitude and Dopplcr modulations that exist in this case. 

Normally, the harmonic at twice the orbital frequency dominates the SNR. However, when the dominant harmonic 
fails to reach the LISA band the higher harmonics become important, which transpires for masses greater than 
4 x 1O 7 M0. Even though the second harmonic falls below the lower cut-off / s of the LISA bandwidth, the kth 
harmonic, k > 2, that has power up to a frequency £;-Flso, might cross / s and produce a significant SNR. Of course, 
the fcth harmonic would fall below the LISA sensitivity band for masses which satisfy the equality / s = fc-FLso- Thus, 
higher PN order waveforms, which bring in higher harmonics, are capable of producing a significant SNR, even when 
the RWF fails to produce any. 

Let us examine this in a little more detail starting from the values of mass where the second harmonic dominates and 
the RWF is adequate. Eventually, for larger values of the total mass, the inequality / s > 3Flso becomes true. Then 
the 0.5PN waveform, which contains the first and the third harmonic, terminates before reaching / s and consequently 
the SNR due to the 0.5PN waveform goes to zero. SNRs for different PN waveforms for a binary whose dominant 
harmonic falls below / s and the third harmonic has a small span in the LISA bandwidth is given in the second column 
of TableU Note that for the 5.5(10 6 - 1O 7 )M system, the 1PN waveform has a higher SNR than the 0.5PN one. This 
is due to the absence of the first harmonic and the small span of the third harmonic in the LISA bandwidth. Further, 
the 2.5PN waveform has a slightly larger SNR compared to 2PN. This is due to the absence of the first and second 
harmonic and the small contribution from the third harmonic, all of which contribute interference terms due to their 
2.5PN corrections. However, this increase is marginal, and is not generic. We have explicitly checked by choosing 
different angles that there can be a small decrease also. The detailed results for LISA are summarised in FigQJ We 
see that for masses for which the 1PN waveform fails to reach the LISA bandwidth, the higher PN order amplitudes 
are capable of producing SNRs as high as 100! Thus, the use of the FWF will enable LISA to make observations of 
SMBHs in the astrophysically interesting mass-regime, which would not be possible had one used only the standard 
RWF. 

Using the expression for Flso, it is simple to argue that the mass reach for the 2.5PN FWF, which has the seventh 
harmonic of the orbital frequency, is 7/2 times the RWF (around 1.5 x 10 8 M Q ). The above ratio, of course, depends 
on the assumption that the Schwarzschild (test particle case) LSO frequency will not be very different from the LSO 
frequency in the comparable mass case. 

We conclude with a discussion of a minor, but clear, feature seen in Fig. [5] for LISA, but not present for the 
ground-based detectors, concerning the relative values of the SNR obtained using the RWF and the FWF. For most 
of the mass range probed the RWF overestimates the SNR relative to the FWF; however, the figure clearly shows an 
anomaly for masses around ~2x 10 6 M Q . To understand this, wc first note that the dips in the two curves in Fig [51 
arc due to the bump in the LISA noise-curve Q just above 10 -3 Hz. This bump is due to the domination of white- 
dwarf confusion noise over instrumental noise and lies just below the most sensitive frequency region (~3x 10~ 3 Hz - 
2 x 10~ 2 Hz) of the LISA band. Below 3 x 10 -3 Hz, the noise increases sharply till one reaches the bump. For binaries 
of mass greater than 1.5 x 10 6 A:/q, the frequency span of the dominant harmonic ends just around the bump and the 
sensitive region of the LISA band is beyond the span of this harmonic. However, higher harmonics incorporated in 
the FWF are able to reach the sensitive part of the noise curve. This leads to higher SNR for the FWF relative to the 
RWF. This reversal of trend continues up to masses 4 x 10 4 A/ Q . Above this mass, the frequency span of the seventh 
harmonic ends before the sensitive region of the LISA band and the general trend is restored. 

For still higher mass values, the SNRs due to the RWF and the FWF both increase until the second harmonic fails 
to reach the LISA band. This is due to the overall scaling of the waveform with the total mass. At such high values, 
it is able to compensate both for the decreasing frequency span and the higher noise of the detector in this frequency 
range. 

E. Effect of higher harmonics in the equal mass case 

In contrast to asymmetric systems discussed so far, for systems of equal mass all odd harmonics are absent. 
Consequently, for symmetric systems the mass-reach of the 2.5PN FWF will be only 3 times the mass-reach of the 
RWF. Further, from the right panel of Fig. [U it is clear that the 0.5PN and the OPN, or RWF, are identical, as 
are the 1PN and 1.5PN waveforms. Thus the decrease in SNR for the higher PN order waveforms with increasing 
total mass is more pronounced than in the unequal-mass case. We also note that for masses for which the second 
harmonic fails to reach the detector bandwidth, the 2PN waveform has a lower SNR than the 2.5PN waveform. This 
can be explained by noting that for these masses only the fourth and sixth harmonics enter the LISA bandwidth. 
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FIG. 5: Luminosity Distance (in Mpc) versus total mass for a fixed SNR of 10. The systems have mass-ratio of 0.1. The 
distance reach can be as large as 500 Gpc for systems where the second harmonic enters the LISA bandwidth. Systems 
undetectable by the RWF (of mass around 1O 8 M0) can be detected by the FWF at distances up to 10 Gpc. The location and 
orientation of the sources are the same as in the earlier figures. 



The 2PN waveform contains the leading term of the fourth harmonic at 1PN and its 2PN correction, which interfere 
destructively. However, inclusion of the 2.5PN amplitude correction leads to a constructive interference term between 
the 2.5PN correction and the 1PN term which is responsible for increasing the SNR for the 2.5PN waveform. 

It is interesting to note that the computation of the 3PN GW polarization which will introduce an harmonic at 8^f 
will be quantitatively more significant for the equal mass case as the mass reach will be better by 33% relative to the 
2.5PN FWF as opposed to the unequal mass case where it is only 14%! This provides one motivation for work in 
progress towards the computation of the 3PN accurate GW polarizations [29[ . 



F. Variation with mass ratio 



Since the mass reach depends only on the total mass, the trends remain the same for different values of mass ratios. 
Fig. [2] compares the variation of SNRs with mass for mass ratios of 0.1 and 0.01. If the SNR is dominated by the 
second harmonic, the SNR is smaller for more asymmetric systems by an overall factor of v, where v = m\ m^lrr? . 
However, once the second harmonic fails to reach the sensitive bandwidth of LISA, the more asymmetric systems have 
a dominant contribution from the odd harmonics which scale by a further factor of \J\ — Av, which is larger for more 
asymmetric systems. Thus the decrease in SNR for the FWF with an increase in the total mass is less steep for more 
asymmetric systems. 



G. Distance reach with the 2.5PN FWF 



Next, we compare the distance-reach of the RWF and the 2.5PN FWF. The results are shown graphically in Fig. [5] 
and are similar in appearance to the mass-reach plot. The mass-reach of the RWF is ~ 4 x 10 7 Ajf Q . For a system of 
total mass 5 x 10 7 M Q , the plot shows that LISA can detect such binaries with an SNR of 10 at a luminosity distance of 
100 Gpc [z ~ 15). SMBHs of total mass ~ 10 8 M Q , not even observable using RWF templates, have a distance-reach 
as high as 10 Gpc (z ~ 1.5) with an SNR of 10. 

Proposals to extend the frequency band-width of LISA up to 10 -5 Hz have been discussed. In that case, the FWF 
can increase the mass-reach of LISA to even around 10 9 M Q . More specifically, LISA can then observe a 10 9 M Q 
system with an SNR of about 30 at 3 Gpc, if it uses templates based on the 2.5PN FWF for data-analysis. 



H. Sensitivity of SNR to source location and orientation 



All the results for SNR using the amplitude-corrected waveforms quoted earlier in this paper have been for a fixed 
choice of location and orientation of the source [defined by the angles(0s, 4>s, 4>l)) with respect to the barycentre 
coordinate system. To conclude our present analysis, in this section we look into the variation in the value of SNR 



0.08 

0.05 I 1 ; 1 1 




SNR SNR due to FWF 

FIG. 6: Distribution of SNR with sources randomly located and oriented in the sky. The left figure plots SNRs due to both 
RWF and FWF for a binary of mass (10 5 - 10 6 )A/ Q . For this mass, the most probable SNR for the FWF is lower than the 
most probable SNR for the RWF, like the trend shown in Table |U The right figure compares the SNRs due to the FWF for 
binaries of mass (10 5 - 10 6 )M o and 5.5(10 6 - 10 7 )M o 



for sources at various locations in the sky and various orientations. To this end, we consider a collection of sources 
randomly oriented in the sky and study the probability distribution of their SNRs. The results of our simulations 
(consisting of 8000 random realisations of the angles involved) arc shown in Fig. [6] From the left panel of Fig. [6] we 
see that the most probable SNR due to the FWF for a (10 5 , 10 6 )A/q binary is less than the most probable SNR due 
to the RWF, indicating that this trend is independent of the source location and orientation. In the right panel we 
see that a binary of mass 2 x (10 6 , 1O 7 )M , which is undetectable by the RWF, can be observed by the FWF with a 
most-probable SNR of around 220. 

V. SUMMARY 

The implications of amplitude corrected 2.5PN full waveforms (FWF) for the construction of detection templates 
for LISA are investigated in detail. With the FWF, LISA can observe sources which are favoured by astronomical 
observations, but not observable with restricted waveforms (RWF). This includes binaries in the mass range 10 8 — 
10 9 A/q, depending on whether the lower cut-off for LISA is chosen to be at 10 -4 Hz or 10~ 5 Hz. With an SNR of 
10, these systems can be observed up to a redshift of about 1.5. The computation of the 3PN polarization, which will 
introduce an harmonic at 8^ (i.e. four times the dominant harmonic), in addition to the existing harmonics, could 
enhance the mass reach for equal mass binaries by 33% and unequal mass binaries by 14.3%. 

The implication of the FWF for parameter estimation will be far more important than the extension of LISA's 
mass-reach reported here. From the work of Van Den Broeck and Sengupta in the context of ground-based detectors 
[24| it is already clear that most parameters will be estimated with errors ~ ten times smaller as compared to RWF. 
This raises the interesting possibility that binary SMBH coalescences might be located on the sky with accuracies 
good enough for optical observations to identify the galaxy cluster and measure its red-shift. Needless to say that this 
improved estimation of source properties will have important consequences in shedding light on the dark energy, better 
understanding of SMBH formation and evolution, structure formation, etc., and is currently under investigation. 

In this work we have confined ourselves to only non-spinning black-holes ignoring the effect of spin-orbit coupling at 
1.5PN [3(j and 2.5PN [3l[ and spin-spin effect at 2PN order [32]. The effect of spin is expected to be astrophysically 
significant and it is important to revisit the present analysis including spin in the future. Though partial results for 
GW polarisations including spin do exist, a more exhaustive exercise would be necessary before the FWF required 
for this work is available. The problem will also be more complicated due to modulations arising from spin-orbit and 
spin-spin couplings which would need to be addressed. 

In this work we also restricted to the inspiral phase and used a physical picture of the LSO that is based on the 
test-particle limit. For comparable masses, the notion of LSO is not as sharp, or unique, and hence our results are 
probably idealized limits of the real situation. Numerical relativity [HI, 1341 13a ] is maturing over the past couple of 
years and could soon provide waveforms for late inspiral and merger. It should then be possible to compare the 
results of such numerical templates with those studied in this paper to provide a better understanding of how higher 
harmonics facilitate the mass reach of our detectors. 
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